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PACS 87.16.D — Membranes, bilayers, and vesicles 
PACS 83 . 50 . Ha - Flow in channels 

PACS 82.70.Uv - Surfactants, micellar solutions, vesicles, lamellae, amphiphilic systems 

Abstract. - The dynamics of fluid vesicles is studied under flow in microchannels, in which the 
width varies periodically along the channel. Three types of flow instabilities of prolate vesicles are 
found. For small quasi-spherical vesicles - compared to the average channel width - perturbation 
theory predicts a transition from a state with orientational oscillations of a fixed prolate shape to 
a state with shape oscillations of symmetrical ellipsoidal or bullet-like shapes with increasing flow 
velocity. Experimentally, such orientational oscillations are observed during the slow migration 
of a vesicle towards the centerline of the channel. For larger vesicles, mesoscale hydrodynamics 
simulations and experiments show similar symmetric shape oscillation at reduced volumes V^* ^ 
0.9. However, for non-spherical vesicles with y* < 0.9, shapes are found with two symmetric or a 
single asymmetric tail. 



Introduction. — Soft deformable objects such as liq- 
uid droplets, vesicles, and cells show a complex behavior 
under flow. For example, in simple shear flow, fluid vesi- 
cles exhibit tank-treading, tumbling, and swinging (also 
called vacillating-breathing, or trembling) motions, de- 
pending on parameters such as shear rate, viscosity con- 
trast, and internal volume [T]-[9]. Understanding the flow 
behavior of lipid vesicles and red blood cells (RBCs) is 
not only an interesting problem of the hydrodynamics of 
deformable, thermally fluctuating membranes, but is also 
important for medical applications. In microcirculation, 
the deformation of RBCs reduces the flow resistance of 
microvessels. In diseases such as sickle cell anemia, RBCs 
have reduced deformability and often block microvascular 
flow [To]. Lipid vesicles are considered as a simple model 
of RBCs and also have applications as drug-delivery sys- 
tems. 

The recent development of microfluidic techniques [11 
allows the investigation and manipulation of individual 
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cells and vesicles, e.g. the measurement of the dynamic 
pressure drop for single-cell deformation [12], the separa- 
tion of RBCs from the suspending plasma [13l|14], and 
the control of oxygen concentration to investigate sickle 
cells [To]. There are many potential applications such as 
blood diagnosis on chip. Under steady flows in homo- 
geneous glass capillaries and rectangular microchannels, 
vesicles [l5] and RBCs [TilfTBHTS] deform into a bullet 
(with a flattened rear end) or parachute (with an inside 
bulge at the rear end) shape. Compared to these steady- 
flow conditions, the vesicle dynamics in time-dependent 
flow is much less explored. Only very recently, phenom- 
ena like the wrinkling of vesicles after inversion of an elon- 
gational flow [19, 2^] or shape oscillation of RBCs [2T1I22] 
and fluid vesicles [23] under oscillatory shear flow have 
been discovered. 

In this paper, we propose a structured microchannel 
system to study vesicle dynamics. The width of a mi- 
crochannel is spatially modulated along the channel, so 
that a flowing vesicle is exposed to an oscillatory elonga- 
tional flow. We have fabricated such microchannels and 
observe the time-dependent vesicle deformation via opti- 
cal microscopy. In parallel, we use perturbation theory 
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for quasi-spherical vesicles and mesoscale-hydrodynamics 
simulations for non-spherical vesicles to predict several 
flow instabilities. In particular, we predict for wider chan- 
nels (compared to the vesicle size) a transition from shape 
oscillations of symmetrical bullet-like shapes to orienta- 
tional oscillations with decreasing flow velocity; for nar- 
rower channels, we observe stable shapes with two sym- 
metric or with a single asymmetric tail. Our results 
demonstrate that structured channels are well suited to 
investigate the dynamical behavior of vesicles. 

Experimental Systems and Methods. — Lipid 
vesicles, consisting of l,2-Dioleoyl-sn-Glycero-3- 
Phosphocholine (DOPC, Avanti Polar Lipids) with 
O.lmol% fluorescently labeled T-Red DHPE (Texas Red 
l,2-dihexadecanoyl-sn-glycero-3-phosphoethanolamine, 
triethylammonium salt, Invitrogen, USA), were prepared 
in aqueous sucrose solution (200mM) by electroformation, 
as described elsewhere [24]. Using this method, many 
giant unilamellar vesicles with diameter larger than 10/im 
were produced. Prior to the experiments, the vesicles were 
introduced into an slightly hypertonic aqueous glucose 
or sucrose solution. This procedure enables control of 
the volume-to-surface ratio by inflation and deflation of 
the vesicles. The vesicles in sucrose solution are perfectly 
density matched. The density difference between sucrose 
and glucose solutions at 200mM is Ap/p = 0.012, which 
implies a small buoyancy force. 

Microchannels were fabricated by soft lithography [301. 
The structure of an AutoCAD designed mask was 
transferred by illumination to a photoresist (SU-8, mi- 
croresist, Berlin), which serves as a master to cast 
poly(dimethylsiloxane) (PDMS) replicas. After treatment 
with an 02-plasma, the PDMS mold was bonded onto a 
thin cover slip. The assembly was connected to a reservoir 
containing the vesicle solution (inlet) and a water reser- 
voir (outlet) to adjust the hydrostatic pressure difference, 
which drives the fluid through the structured microchan- 
nel (see Fig. [T|. All channels have a height Lz = 50/im 
and a periodic length = 100/im. 

Due to the lift force from a wall [25H27) and the mi- 
gration in a parabolic Poiseuille flow [SS' 29 , vesicles are 
aligned close to the center of the channel. In the glu- 
cose solutions, the buoyancy force leads to a displace- 
ment of the vesicles from the center of the channel in 
the z-direction (see Fig. [1]). The balance of the buoy- 
ancy force and the lift force [25l[27] implies a distance 
I {r]jRo/{Apg)y^'^ between membrane and wall, where 
7 is an effective shear rate, Rq a mean vesicle radius, rj 
the fluid viscosity, and g the gravitational acceleration. 
For Rq = 10/im, 7 = ls~^ (corresponding to a mean flow 
velocity of about 20/im/s) and the viscosity of water, we 
obtain the order-of-magnitude estimate of = 10/im, com- 
parable to the channel size L^. The displacement from the 
center implies an asymmetric deformation of the vesicle 
shape in glucose solution under flow. 




Fig. 1: (Color online) Side view of the channel setup with 
the reservoirs. The inset shows a top view of the modulated 
channel. The length of the modulation is Lx = 100/im. 

Theory for quasi-spherical vesicles in modulated 
channels. — We consider vesicles which have the same 
viscosity rj of the outer and inner fluids. Vesicles have 
constant volume V and constant surface area so that 
the reduced volume V* and the excess area As are defined 
by = {Rv/Rsf = (l + As/4^)-^/2 and As = S/R^- 
47r, where Ry = {SV/Air)^^^ and Rs = (6'/47r)i/^ 

First, we derive an analytical description of the flow be- 
havior of small quasi-spherical vesicles with As ^ 1, based 
on the theory for linear flows [21I5H71I2Q] . The microchan- 
nel has a constant height and a periodically modulated 
width with walls at :^hy{x)^ where 

hy{x) = ^ {1 + ay cos [^)}. (1) 

Here, Lx is the periodicity length along the channel (see 
Fig.[T]). The flow field v(r) for weakly modulated channels 
{ay <^ Lx/Ly) is obtained from lubrication theory, 

Vx = {9v^Ly/Shy){l-{y/hyf}{l-{2z/Lzf}, 
Vy = -7rVxyLyaySm{27rx/Lx)/hyLx, (2) 
Vz = 0, 

where Vm is the mean flow velocity. Furthermore, we con- 
sider small vesicles with Ry <C Ly^Lz^ and neglect back- 
flow effects due to the channel wall. 

The vesicle shape is expanded in spherical coordinates 
as r = Ry{l + ^ 'U|,m^i,m), with the spherical harmon- 
ics Yi^mi^^^) and the polar axis in the z direction. We 
consider a vesicle shape which is mirror symmetric in the 
z direction, which implies, e.^., ii2,±i = '^3,0 = '^3, ±2 = 0. 
The curvature energy of the membrane with bending rigid- 
ity tv and reduced surface tension a is given by 

F = {K/2)Y,Ei\ui^m\^ + k{U + 2a/3)4'\ (3) 

where El = (/ + 2)(/- l){/(/ + l) + cr}. Here, the third order 

term 1^2^^ = (y^S/tt/T) (1^2,0^ — 6ii2,o|'^2,2p) for / = 2 is 
necessary to obtain a prolate shape as thermal-equilibrium 
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state [71I3T] . With the Stokes approximation and the Lamb 
solution for the flow field, the dynamics of the amplitudes 
ui^rn of a vesicle moving along the center line {y = z = 0) of 
the channel with velocity v(xg,0,0) = (9'UmL^/8/iy, 0, 0) 
is governed by 



dt 



(4) 



with reduced bending rigidity K 
Ly/{vmr) and = /(/ + l)/(2/ + l)(2/2 



I^Ly/{r]Ry^Vn 



.) 



21 - 1). Here, 



r = r]Ry / Hi is the characteristic relaxation time of a vesi- 
cle. Two kinds of forces determine the membrane defor- 
mation in Eq. (|4]), the curvature force fi^rn = Eiui^rn — 
(12 + 2a/3)du2 /dui^rn^ and the force si^rn{xQ) due to 
the elongational fiow. The force s^^^(xg) is given by 
si,m{xQ) = {ai^rn ^ Bibi^rn)Ly / RyVm [lEllin] • Here, Bi = 
l/(2/2+2/-l),and 

<^Lm and bi rn are the spherical harmon- 
ics expansion of Vr{r) and Rydvriv) / dr at |r — tqI = Rw^ 
respectively, where Vr(v) is the radial components of the 
fiow field v(r) in the absence of the vesicle. The tension 
cr is determined by the area constraint dAs/dt = with 
As = Ez,m(^ + - l)\ui^rn\V2 - 2uf^ We cousidcr a 
vesicle with a small excess area As = 0.1 (corresponding 
to V* = 0.988), where it is sufficient to take into account 
/ = 2, 3 modes, and a channel with Ly = = 20Ry and 
various Lx and ay. 

In fast fiow, where fiow forces dominate and the reduced 
bending rigidity is small, a^* <C 1, the analysis of Eq. (|4]) 
shows that the shape is mirror symmetric with respect 
to both xz and xy planes, and that the vesicle lengths 
lx and ly oscillate (see Fig. [2]). The fiow elongates the 
vesicle in the y and x direction for —Lx/2 < x < and 
< X < Lx/2, respectively. In the regime z^* <C 1, the 
extremal elongations are determined by the balance of the 
fiow forces, si^rnixc)^ and the area constraint. The forces 
1^52,2(^0) I and 1 53^3 (xg) I have maxima at xq = ±0.38^^; 
and xq = ^Lxl2 for = 0.5, respectively. The resulting 
positions of extremal elongations are found to be close to 
xq = and xq = Lxl2 (see Fig. Efb)). 

In contrast, in slow fiow with ^ \, the vesicle is 
predicted not to change its shape, but to display a periodic 
oscillation of its orientation. The tilt angle d oscillates 
around 7r/4 (or — 7r/4 depending on initial positions), see 
snapshot in Fig. [2fa). Thus, a symmetry breaking occurs 
with decreasing fiow velocity. 

Both types of vesicle motions are limit cycles. A sharp 
but continuous transition between these two motions oc- 
curs at a critical reduced bending rigidity a^*, where the 
tilt angle d vanishes (see Fig. EJa)). The critical value 
increases with increasing corrugation ay of the chan- 
nel, but is almost independent of Lx for Lx/Ly > 2 (see 
Fig.[2fc)). The symmetric vesicle deformation reduces the 
disturbance of the original fiow but increases the free en- 
ergy of a vesicle. At small or large the former or latter 
contribution dominates, respectively. 
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Fig. 2: (Color online) Dynamics of a quasi-spherical vesi- 
cle at = 0.988. (a) Dependence of the tilt angle Q 
of a vesicle on ac*, for Lx/Ly — 4 and ay — 0.5. Here, 
(sin^(2^)) = ((Im(i^22)/ 1^^22 1)^) describes the deviation from 
symmetric shape (with respect to the xz plane). The insets 
show sliced snapshots of vesicles in the xy plane for ^c* = 0.01 
and AC* = 50. Solid and dashed lines indicate shapes of ex- 
tremal elongation or tilt, (b) Maximum vesicle lengths in the 
x,y,z directions as a function of the center-of-mass position 
xg, for K* = 0.01 and ay = 0.5. Solid and dashed lines rep- 
resent results for Lx/Ly — 4 and Lx/Ly — 16, respectively, 
(c) Critical reduced bending rigidity Kjq as a function of the 
corrugation amplitudes ay for channel geometry Lx/Ly — 2, 
4, 8, and 16. The inset shows H\jQ as a function of Lx/Ly for 



Mesoscale Hydrodynamics Simulations of Vesi- 
cles under Flow. — We employ mesoscale hydrodynam- 
ics simulations to study the flow behavior of vesicles in nar- 
row channels. A dynamically-triangulated surface model 
for the membrane combined with a particle-based 

mesoscale simulation technique — multi-particle collision 
dynamics (MFC) j33H35] — for the embedding fluid. A 
detailed description of this approach to model vesicles dy- 
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Fig. 3: (Color online) Sequence of snapshots (at equal time 
intervals) of vesicles with = 0.96 and Rs/Ly = 0.21 moving 
through a microchannel, experimentally observed by optical 
microscopy (in sucrose solution) for = 0.04 {vm = 38/im/sec, 
Ly — 75/im) (top), and in simulation for k* = 0.08 (bottom). 
For the experimental data, the time interval between images is 
0.97 sec. 



namics under flow can be found in Refs. [4l|T8]. In MFC, 
the fluid is described by point-like particles of mass 
mg. After free streaming of every particle for a time 
step At, the fluid particles collide in cubic boxes of lat- 
tice constant a = L^^/IG, which is also the mean length 
between membrane vertices. We use the fluid viscosity 
T] = bb{)\/m^kBT / (with number density n = 100a~^ 
and the time step At = O.Olay^m^T^B^), corresponding 
to low Reynolds numbers, to simulate experimental condi- 
tions. We study vesicles with bending rigidity n = 20/cbT 
where k^T is the thermal energy. The flow velocity is 
chosen to be v^t / Ly = 17 {) {Rsj / Ly^ for a channel with 
Ly = Lz = Lx/2, corresponding to a slow vesicle relax- 
ation compared to the passage time through a channel seg- 
ment. For vesicles with Rs/Ly = 0.36 and Rs/Ly = 0.46 
at = 0.9, the reduced bending rigidity is hz* = 0.14 
and = 0.07, respectively. In this regime, the vesicle 
dynamics is not very sensitive to a^*. 




Fig. 4: (Color online) Maximum lengths Ix, ly, and lz of vesi- 
cles in experiments (solid lines) and simulations (dashed lines) 
from the same data as shown in Fig. \3\ 

Vesicles in sawtooth-shaped channels in fast 
flows. — We investigate the dynamics of large prolate 
vesicles in fast flows, corresponding to small a^*, both by 
experiments and simulations. In this case, hydrodynamic 
interactions between the vesicle and wall are not negligi- 
ble. We consider now a periodically-patterned microchan- 
nel with a saw-tooth shape (see Figs. [T]and[3j), for which 
the walls are located at 

hyix) = ^{l + ay[l±^)} (5) 

for —Lx/2 < X < and < x < Lx/2^ respectively. 
Two types of channels are used, which are characterized 
by Ly = 50/im, ay = 0.2 and Ly = 75/im, ay = 0.33 
(as well as Lz = 50/im and L^ = 100/im). The glucose 
and sucrose solutions are used for the experiments with 
the narrower and wider channels, respectively. Typical 
experimental flow velocities are in the range Vm = 10 - 
100/im/s. 

First, we describe the dynamics in the wider channel 
{Ly = 75/im). Experimental and simulation results for 
the shape deformation are shown in Figs. [3] and IH As ex- 
pected, the vesicles change their shapes periodically. For 
^0.9, vesicle shapes vary continuously between an al- 
most ellipsoidal (in the narrow part of the channel) and a 
cone-like (in the wide part of the channel) bullet shape, see 
Fig. [3l This is in contrast to the behavior for unbounded 
steady Poiseuille flow, where perturbation theory predicts 
a coexistence between bullet- and parachute-like shapes at 
the center line [29^ . The amplitudes and phases of the de- 
formations in the x- and ^/-directions (see Fig.H]) agree well 
between simulations and experiments, while very little de- 
formation is seen in the z-direction. The small difference 
of amplitudes between the experiment and simulation is 
probably due to the somewhat different values of the re- 
duced volume and radius of the vesicle, since an accurate 
experimental estimation is difficult (because the length in 
the z direction could not be measured). The dynamics 
also agrees qualitatively with the theoretical predictions 
for quasi-spherical vesicles, compare Fig. [21 although the 
parameters are outside the range of validity of the theo- 
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Fig. 5: (Color online) Normalized oscillation amplitudes Al/{1) 
in the x (red) and y (blue) directions. Open symbols represent 
simulation data for Rs/Ly = 0.36 (o, □) and Rs/Ly = 0.46 
(A, o) at tty = 0.2, respectively. Closed symbols represent 
experimental data for Ly = 50/im with ay = 0.2 in glucose 
solution (•, ■) and Ly = 75/im with ay = 0.33 in sucrose solu- 
tion (a, ♦), respectively. The solid lines are guide to the eye. 
The lower part shows snapshots of simulations (si, s2, s3, s4) 
and experiments (el, e2, e3, e4) for Ly = 50/xm. All snapshots 
are displayed with the same length scale. The pictures el and 
e2 show the same vesicle before and after the transformation 
from ellipsoid to tailed shape, respectively. All images (el-e4) 
are vesicles in glucose solution. 

retical approximations. 

At reduced volumes < 0.9, we find that the vesicles 
deform into novel shapes, while they are still ellipsoidal in 
the absence of flow. At V* ^ 0.85, the simulated vesicles 
develop two tails in the xy plane (see Fig.[5fs3)). This tail 
position is very stable; the tails quickly return to the xy 
plane, even from an initial state with the tails in xz plane 
(obtained by 7r/2 rotation). At V* 0.80, the symmetric 
tails become unstable; they are replaced by a stable shape 
with a single asymmetric tail (see Fig.[5fs2)). In this case, 
the configuration with a tail in the xz plane is also stable 
within the accessible simulation time (see Fig. [5jsl)). A 
single-tail shape is also observed in our experiments with 
vesicles in glucose solution (see Fig. [He2)); here, the off- 
center motion induced by the buoyancy force implies that 




Ol I I I I I I I I I I I I I L 



° ^Xq/L, 2 3 

Fig. 6: (Color online) Oriental oscillation of a vesicle (in su- 
crose solution) in a wide structured channel with Ly — 75/jjm. 

(a) Micrograph of an oscillating vesicle with Ry = 24/im, 
1' = 34 /xm/s, corresponding to — 0.0132. The image is 
a superposition of video frames taken at time intervals of Is. 

(b) Orientational angle as a function of the center-of-mass po- 
sition xg- The (red) circles (o) represent the same data as in 
(a). The vesicle oscillates between the extrema of ^min = O.OTtt 
and ^max = 0.37r. The (blue) triangles represent another vesi- 
cle (with Rv = 18/im, v = ib55/im/s). Here, we have inverted 
the direction of the fluid flow; symbols indicate angles before 
(A) and after (a) flow inversion. 

the tail is found in the z direction. Recently, Abkarian et 
al. [H] found a similar asymmetric tail for RBCs in very 
fast flows in homogeneous capillaries (with flow velocity 
'Um = 10 — 35 cm/s and capillary diameter 10/im). Such 
a long tail is much more difficult to form in RBCs due 
to the shear resistance of the spectrin network, and may 
even require a local separation of the lipid bilayer from the 
spectrin network. Thus, the formation of an asymmetric 
long tail is a generic feature in capillary flows. Structured 
channels promote tail formation already at smaller flow 
velocities. The low contrast and the blurry contour in the 
rear part of the vesicles in the experimental micrographs 
of Fig. [5] (e3, e4) indicates that vesicles exhibit enhanced 
shape fluctuations at the rear part, due to a locally re- 
duced membrane tension; the low contrast regions in the 
shapes of Fig. [5] (el, e2) is probably due to a tilt of the 
vesicle axis and a vesicle asymmetry induced by off-center 
motion with a buoyancy force. 

Our experimental and simulation results for the oscilla- 
tion amplitude A/, displayed in top part of Fig. [5l show 
that the difference between maximum and minimum of 
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Ix or ly is quite insensitive to the reduced volume 
and the vesicle size Rs. In the experiments, V* is esti- 
mated from the prolate ellipsoid shape with the lengths 
{Ix) and {ly) of the long and the two short axes, respec- 
tively. The finite amplitudes for V* 1 are caused by 
errors in the estimation of F*. The velocity v^q^ and av- 
erage size Im = {{Ix) + (^a?))/2 of liposomes are varied in 
experiments in the ranges v^es = 14 - 45/im/s and Im = 20 
- 40/im, respectively. The amplitude Aly of the vesicle 
with asymmetric tail in the xy plane (s2) is smaller than 
in xz plane (si), since its tail is still at the wide part of 
the channel and elongated when the main body is at its 
narrowest part. 

For vesicles, which fiow in the center of the channel, sim- 
ilar symmetric shape oscillations are observed in the nar- 
row and wide channels (in agreement with our predictions 
for fast fiows from perturbation theory). However, when 
the vesicle is slightly displaced from the center {y = 0) in 
the wider channel it displays orient ational oscillations, as 
shown in Fig.[6l Here, the asymmetry of the flow field acts 
to magnify the orient ational oscillations, which are accom- 
panied by strong shape oscillations between slipper and 
ellipsoidal shapes. Thus, a mixed state of orientational 
and shape oscillations appears for /^* < 1; it appears to be 
long-lived due to the small lift force for wider channels. 

Vesicles in sawtooth- shaped channels in slow 
flows. — Finally, we want to briefly discuss the dynam- 
ics of large vesicles at slow flows, with a^* ^ 1. If the 
vesicle size is larger than the smallest channel width L^^", 
vesicles can not flow through the channel without a shape 
deformation. Thus, with decreasing flow velocity, the dy- 
namics changes from shape oscillations to a trapped state 
[36] (instead of the orientational oscillation). The tran- 
sition velocity depends on Rs/L^^^^ F*, and a^*. A mi- 
crochannel with a strong variation of channel width (for 
instance the channel in Fig. 6(a) of Ref. 114^), is suitable 
for a study of this trapped-escape transition. 

Summary. — We have studied the dynamics of vesi- 
cles in structured microchannels, both theoretically and 
experimentally. For large reduced volumes, the vesicles 
periodically change their shape (large flow rates) or their 
orientation (small flow rates), depending also on the chan- 
nel geometry. For smaller reduced volumes, we flnd a novel 
shape with a long asymmetric tail. The good agreement 
of theoretical and experimental results shows that flow of 
vesicle suspensions in complex flow geometries can now be 
understood quantitatively. 

* * * 
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